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CO  ENERGY  TRANSFER  PROCESSES  IN  THE 
MESOSPHERE  AND  THERMOSPHERE 


The  original  purpose  of  this  review  was  to  update  the  CO  chemical  kinetics  module 
in  the  Strategic  High  Altitude  Radiance  Code  (SHARC)  (Ref.  1).  The  code  contained  circa 
1974  rate  constants  for  the  collisional  relaxation  of  vibrationally  excited  CO.  The  source  of 
these  rate  constants  was  a  review  article  by  Taylor,  entitled  Energy  Transfer  Processes  in 
the  Stratosphere  (Ref.  2). 

Since  1974,  new  experimental  techniques  have  enabled  researchers  to  study  the 
temperature  dependence  of  specific  vibrational  energy  relaxation  processes.  In  particular, 
laser-induced  fluorescence  (LIF)  studies  have  filled  the  gap  between  room  temperature 
measurements  and  high  temperature  shock-tube  studies.  The  data  acquired  for  CO  has 
improved  in  both  quantity  and  reliability  over  the  past  few  decades. 

For  the  purpose  of  this  report  and  for  updating  SHARC,  CO  relaxation  processes 
relevant  to  the  thermosphere  and  mesosphere  have  been  considered.  The  altitude  range 
covered  by  the  model  atmosphere  profiles  in  SHARC  (Version  2.0,  1991)  is  50  to  300  km. 
Over  this  altitude  span,  the  density  of  molecules  decreases  dramatically,  so  the  number  of 
collisions  also  decreases.  For  vibrationally  excited  species  that  are  potential  infrared 
radiators,  spontaneous  emission  processes  become  competitive  with  nonradiative 
collisional  deactivation;  that  is,  the  system  begins  to  deviate  from  local-thermodynamic 
equilibrium  (LTF).  Generally,  if  there  are  less  than  a  million  collisions  in  a  radiative 
lifetime  1/A,  where  A  is  the  Einstein  coefficient  for  spontaneous  emission,  then  the  system 
is  non-LTE  (Ref.  3).  In  order  to  fully  characterize  the  atmospheric,  non-LTE  processes  of 
CO,  the  Einstein  A  coefficients  for  each  relevant  vibrational  level  have  also  been  reviewed. 

The  available  literature  (through  1990)  was  evaluated  for  each  of  the  following 
collisional  deexcitation  pathways,  where  v  =1  or  2: 


CO(1)  +  N2(0)  =  CO(0)  +  N2(1)  (1) 

CO(  1 )  +  0^0)  =  CO(0)  +  0^1)  (2) 

CO(v)  +  O  =  CO(v  -  1  or  v  -  2)  +  O  (3) 

CO(v)  +  M  =  CO(v  -  1  or  v  -  2)  +  M,  M  =  N2  or  02  (4) 


In  the  case  of  02  and  N2  as  collision  partners,  both  V-V  and  V-T  relaxation  mechanisms 
were  considered.  The  relative  spacing  of  the  vibrational  energy  levels  of  CO,  02.  and  N2 
is  illustrated  in  Figure  1. 
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Figure  1.  Relaxation  mechanisms  of  vibrationally  excited  CO. 


Ambient  temperatures  between  approximately  190  and  1300  K  are  common  10  me 
altitude  range  from  50  to  3(X)  km.  The  lowest  temperatures  tire  likely  to  be  found  in  the  85- 
90  km  altitude  layer  at  arctic  latitudes,  while  the  highest  temperatures  that  occur  at  the  top  of 
the  atmosphere  (300  km  in  SHARC)  are  produced  during  episodes  of  high  solar  activity. 
A  common  functional  form  for  the  temperature  dependence  of  the  rate  constants  for 
collisional  deexcitation  (and  excitation)  reactions  is  desirable  for  modeling  purposes.  The 
following  form,  used  by  SHARC  and  adopted  in  this  review,  is  based  on  the  theory  of 
Landau  and  Teller  (Ref.  4): 

k  =  Aj  exp  (-Bj/T1^  -  Q/T)  (5) 

The  T'03  dependence  arises  from  the  expression  for  the  most  probable  velocity  for  an 
effective  energy  transfer  collision  (Refs.  4,  5).  In  the  present  analysis  of  experimental  rate 
constant  data,  plots  of  In  k  vs.  T-1/-3  were  used  to  establish  linear  trends  that  can  be 
extrapolated  or  interpolated  to  temperatures  of  atmospheric  relevance.  The  -Cj  /T  factor  in 
Eq.  (5)  was  added  to  the  exponential  term  to  describe  processes  that  have  an  energy  defect. 
If  the  energy  transfer  process  is  endothermic.  C,  =  AE:  otherwise,  Cj  =  0.  The  collisional 
relaxation  data  are  compiled  in  Appendices  A-D.  The  original  references  are  listed  with  the 
appendices. 

The  spontaneous  emission  processes  that  can  occur  from  the  CO  v=l  or  v=2 
vihnuinnal  levels  are 


CO(  1 )  — >  CO(0)  +  hv 

(6) 

CO(2)  — >  COO)  +  hv 

(7) 

CO(2)  CO(0)  +  hv 

(8) 

The  unimolecular  rate  constants  for  these  decay  processes  can  be  calculated  from 
the  appropriate  rovibrational  band  strength  or  the  transition  dipole  moment.  In  order  to 
standardize  all  of  these  rates,  the  most  recent  sum-of-lines  band  strengths  from  the 
Atmospheric  Trace  Molecule  Spectroscopy  (ATMOS)  (Ref.  b)  high  resolution  data  base 
have  been  used  to  calculate  the  Einstein  A  coefficients  for  the  processes  described  by  Eqs. 
(6)  (8).  These  values,  as  well  as  references  to  the  spectroscopic  literature,  are  given  in 
Appendix  E. 

Table  I  summarizes  the  recommended  changes  to  the  rate  constants  for  both  the 
radiative  and  nonradiative  processes.  CO  radiance  calculations  were  performed  with 
SHARC  using  both  the  old  (1974)  collisional  relaxation  rates  and  the  revised  (1990) 
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(d)  Rate  units  lire  cm!  molecule-1 

(e)  Rate  units  are  s'1. 


collisional  relaxation  rates.  The  small  changes  made  to  the  Einstein  A  coefficients  (see 
Table  1)  were  not  included  in  the  calculation.  An  equatorial  atmosphere  that  represented 
relatively  high  solar  and  geomagnetic  activity  was  constructed  using  the  Mass  Spectrometer 
and  Incoherent  Scatter  (MS1S)  1983  model  (Ref.  7)  and  the  SHARC  1976  Standard 
Atmosphere  (Ref.  If  This  particular  atmosphere,  which  features  high  O  atom  densities  at 
altitudes  >  86  km,  was  chosen  for  its  sensitivity  towards  changes  in  the  CO  kinetics.  A 
listing  of  the  atmosphere  is  provided  in  Appendix  F.  Using  this  model  atmosphere  for  the 
starting  molecular  number  densities,  SHARC  calculated  excited  state  populations  for  CO. 
The  other  species  included  in  the  atmosphere  were  H,  O,  NN,  Og,  ICO,  OH.  NO,  COj, 
and  Ov  The  results  of  the  comparison  are  given  in  Table  2.  The  overall  effect  of  the 
revisions  was  an  increase  in  the  intensity  of  the  ( 1-0)  band  at  the  expense  of  the  '2  1  i  and 
(2-0)  bands.  The  total  band  radiance  for  the  three  bands  was  essentially  constant  Since 
the  (2  ! )  and  (2-0)  are  much  weaker  than  the  ( 1-0)  band,  a  factor  of  100  change  in  the  ho; 
and  overtone  bands  resulted  in  only  a  10%  difference  in  the  ( 1  -0)  band.  The  el  feet  on  the 
CO  (2-0)  emission  band  is  illustrated  in  Figure  2.  The  diagram  compares  the  radiance 
calculated  by  SHARC  using  the  set  of  rate  constants  proposed  by  Taylor  (Ref.  2)  with  that 
calculated  using  the  rate  constants  listed  in  this  review  (Table  1 ). 

The  experimental  data  for  CO  relaxation  processes  have  considerably  improved 
since  1974.  The  updated  temperature  dependences  calculated  in  this  report  should  be 
included  in  radiance  models  used  for  infrared  sensor  designs  covering  the  4.7  pm  and/or 
the  2.3  pm  regions.  Energy  transfer  rates,  including  those  for  CO.  are  also  important 
parameters  in  climate  models  used  to  calculate  the  dissipation  of  sunlight  or  electron  energy 
entering  the  Earth's  atmosphere.  Absorption  and  reemission  of  radiation  by  atmospheric 
molecules  is  an  important  consideration  in  calculating  global  temperature  trends. 


If) 


Table  2.  CO  Band  Radiance  Comparison^Mb) 

A  B 

CO  Vibrational  Band  Radiance  Radiance  Comparison 

(Taylor,  1974)  (Revised  Kinetics)  A/B 

W/sr/cm2  W/sr/em2 


(1-0)  2143.27  cm1 

1.6456e-10 

1.81 29e- 1 0 

0.91 

(2-1)  2116.79  cm-1 

L6378e-1 1 

1.6529e-13 

99.1 

(2-0)  4260.06  cm-1 

7.0440e-13 

7.0969e- 1 5 

99.3 

TOTAL 

1.8164e-09 

1.8146e-09 

1.00 

(a)  See  Appendix  F  for  Model  Atmosphere  Listing. 

(b)  SHARC  2.0  Calculation,  200  km  Limb  View,  0°  Latitude,  0°  Longitude. 
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Figure  2.  Calculation  of  C0(2-G)  radiance  using  SHARC  2.0. 
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APPENDIX  A:  V-V  ENERGY  TRANSFER  FROM  CO  TO  02 


In  the  forward  direction,  the  process  described  by  Eq.  (Al) 

CO(  1 )+  02(0)  =  CO(O)  +  02(  1 )  (A  1 ) 

is  an  exothermic  energy  transfer  process.  The  energy  discrepancy  of  587  cm'1  between  the 
C0(1)  and  02(1)  levels  is  too  large  to  classify  this  process  as  resonant,  even  though  it 
occurs  between  discrete  vibrational  energy  levels  of  two  molecules. 

Taylor  (Ref.  Al)  derived  the  following  temperature  dependence  for  Eq.  (Al)  based 
on  high  temperature  shock-tube  references  (Refs.  A2,  A3): 

k  forward  =  3.50e-10  exp  (-124.0  T'1/3)  (A2) 

Several  more  recent  temperature  dependent  studies  have  been  done.  Doyennette  et 
al.  (Ref.  A4)  and  Gregory  et  al.  (Ref.  A5)  performed  laser-induced  fluorescence  (LIF) 
studies,  while  Miller  and  Millikan  (Ref.  A6)  used  a  steady  state  vibrational  quenching 
technique.  The  more  limited  LIF  room  temperature  studies  of  Green  and  Hancock  (Ref. 
A7)  and  the  chemiluminescence  studies  of  Hancock  and  Smith  (Ref.  A8)  overlap  with  the 
data  of  Refs.  A4-A6.  The  data  from  Refs.  A2-A8  are  compiled  in  Figure  A 1 .  The  Landau- 
Teller  expression  obtained  from  the  slope  and  intercept  of  Figure  Al  is 

k forward  =  1.17e-12  exp(-61.3  T-!/3)  (A3) 

The  appropriate  reverse  rate  equation,  including  the  energy  difference  term,  is 

simply 

krevcrsc  =  117e-12  exp(-61.3  T  -  844.4  T"1  )  (A4) 

There  is  a  significant  difference  between  the  temperature  dependence  described  by 
Taylor  and  the  temperature  obtained  in  this  analysis,  according  to  Figure  Al.  The  CO  to 
02  V-V  rates  determined  experimentally  by  Refs.  A4-A8  are  generally  faster  than  those 
previously  predicted  by  Taylor's  extrapolation  of  high  temperature  shock-tube  data  to  low 
temperatures. 

At  temperatures  below  ~  600  K,  each  of  Refs.  A4-A6  shows  a  steep  curvature  in 
the  Landau-Teller  plot.  Since  the  curvature  is  very  different  in  each  set  of  data,  the  effect  is 
averaged  out  in  the  linear  fit  of  Eq.  (A3).  References  A4  and  A5  contain  a  discussion  of 
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t-1/3  (k-1/3) 
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0.3 


ergy  transfer  from  CO(  1 )  to  O2  (0).  Data  from 
1  (Ref.  A6),  □ ;  Doyennette  et  al.  (Ref.  A4),  • ; 

1  (Ref.  A8)/Green  and  Hancock  (Ref.  Al),  A; 

1 ;  Gregory  et  al.  (Ref.  A5),  A  ;  Sato  et  al.  (Ref. 
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the  importance  of  short-range  intermolecular  forces  in  this  type  of  energy  transfer  process. 
Over  narrower  temperature  ranges,  different  functions  of  T  could  provide  a  better  fit  to  the 
data. 


17 


REFERENCES 


Al.  R.  L.  Taylor,  Can.  J.  Chem.,  52,  1381  (1974). 

A2.  Y.  Sato,  S.  Tsuchiya,  and  K.  Koratini,  J.  Chem.  Phys.,  50,  1911  (1969).  (1000- 
3000  K) 

A3.  R.  E.  Center,  J.  Chem.  Phys.,  58,  5230  (1973).  (1500-4000  K) 

A4.  L.  Doyennette  et  al.,  J.  Chem.  Phys.,  67,  3360(1977).  (625-100  K) 

A5.  E.  Gregory  et  al.,  J.  Chem.  Phys.,  78,  3881  (1983).  (90-300  K) 

A6.  D.  Miller  and  R.  Millikan,  Chem.  Phys.  Lett.,  27,  10  (1974).  (300-140  K) 

A7.  W.  H.  Green  and  J.  K.  Hancock,  J.  Chem.  Phys.,  59,  4326  (1973).  (Room 
Temperature) 

A8.  G.  Hancock  and  I.  W.  M.  Smith,  Appl.  Opt.,  10,  1827  (1971).  (Room 
Temperature) 


19 


APPENDIX  B:  V-V  ENERGY  TRANSFER  FROM  CO  TO  N2 


The  following  process  is  a  resonant  energy  transfer  since  the  CO(l)  level  occurs 
only  187.5  cm'1  below  the  14N2(1)  level: 

CO(l)  +  N2(0)  =  CO(O)  +  N2  (1)  (Bl) 

The  equation  used  by  Taylor  to  describe  the  temperature  dependence  of  the 
endothermic  process  was  (Ref.  Bl) 

^forward  =  6.98e-13  exp  (-25.60T'1/3  -  268.5T'1)  (B2) 

The  rate  equation  for  the  reverse  reaction  has  the  same  coefficients  without  the  AE/T  term. 
Four  references  (B2-B5)  were  cited  by  Taylor.  Of  these,  only  one  room  temperature  study 
by  Zittel  and  Moore  (Ref.  B2)  employed  the  laser-induced  fluorescence  (LEF)  technique. 
The  others  were  high  temperature  shock-tube  measurements.  However,  due  to  this  one 
reliable  room  temperature  measurement,  the  temperature  dependence  that  Taylor  derived 
(Eq.  (B2)J  was  not  appreciably  affected  by  later  measurements  at  temperatures  between  70 
and  700  K.  In  the  present  evaluation,  the  data  from  Refs.  B2-B5  were  combined  with 
newer  data  from  Refs.  B6-B1 1.  In  Figure  Bl,  In  k  +  268.5  T'1  vs.  T'1/3  is  plotted  in  order 
to  directly  fit  for  the  temperature  dependence  of  CO  deactivation  by  N2.  A  distinct  change 
in  slope  occurs  at  T"1/3  =  .08  (i.e.,  T  =  2000  K).  Neither  Taylor's  Eq.  (Bl)  nor  the 
present  fit  described  by  the  following  equation 

^forward  =  9.06e- 1 3  exp  (-26.84T'1/3  -  268.5T'1)  (B3) 

includes  the  measurements  for  T  >  2000  K,  since  temperatures  below  2000  K  are  generally 
used  for  atmospheric  calculations. 

On  an  expanded  T"1/3  scale,  each  set  of  low  temperature  results  deviates  positively 
from  the  linear  fit,  particularly  for  T  <  200  K.  A  discussion  of  the  theoretical  basis  for  the 
observed  temperature  trends  in  the  rate  of  V-V  energy  exchange  between  diatomic 
molecules  can  be  found  in  Ref.  B6.  Maricq  et  al.  describe  a  Morse  intermolecu'ar  potential 
that  can  be  used  to  adequately  describe  the  behavior.  However,  the  low  temperature 
curvature  away  from  Landau-Teller  behavior  appears  to  be  less  than  Uie  error  limits  due  to 
systematic  variations  between  individual  experiments.  For  example,  the  data  of 
Mastrocinque  are  somewhat  inconsistent  with  the  shape  reported  in  the  other  references. 
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Generally,  the  curvature  becomes  well  defined  at  temperatures  that  are  lower  than  those 
relevant  to  mesospheric  processes,  so  we  have  maintained  the  linear  function  in  T'1^. 
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Figure  Bl.  V-V  energy  transfer  from  CO(l)  to  N2(0).  Data  from 
Mastrocinque  et  al.  (Ref.  B8),  ■;  Stephenson  and  Mosburg  (Ref.  B10),+; 
McLaren  and  Appleton  (Ref.  B4),  □;  Maricq  et  al.,  (Ref.  B6),0;  Zittel 
and  Moore  (Ref.  B2),  /Green  and  Hancock  (Ref.  B9),»;  Starr  et  al.  (Ref. 
Bl  !),•;  Sato  et  al.  (Ref.  B3),x;  von  Rosenberg  et  al.  (Ref.  B5),A. 
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APPENDIX  C:  V-T  ENERGY  TRANSFER  FROM  CO  TO  Q(3P) 


At  altitudes  above  approximately  100-200  km,  depending  on  the  atmospheric 
model,  ground  state  atomic  oxygen  becomes  the  dominant  species.  In  some  cases,  such  as 
the  atmosphere  described  in  Appendix  F,  there  is  a  relatively  high  O  atom  number  density, 
even  in  the  mesosphere.  CO  vibrational  and  translational  relaxation  by  ground  state  atomic 
oxygen  can  occur  via  the  following  pathways: 


CO(l)  +  0(3P)  =  0(3p)  +  CO(0) 

(Cl) 

CO(2)  +  0(3P)  =  0(3P)  +  COG) 

(C2) 

CO(2)  +  0(3P)  =  0(3P)  -i-  CO(0) 

(C3) 

The  temperature  dependence  of  the  rate  of  Eq.  (Cl),  derived  by  Taylor  (Ref.  Cl) 
in  1974,  was 


ki  forward  =  9.90e-8  exp  (-118.1T*'/3)  (C4) 

This  expression  was  based  entirely  on  high  temperature  shock-tube  data  (Refs.  C2,  C3). 

More  recent  laser-induced  fluorescence  (LIF)  data  between  265  and  389  K  (Ref. 
C4)  have  been  added  to  extend  the  temperature  range  and  get  an  improved  fit  for  kj.  The 
data  from  Refs.  C2-C4  are  plotted  in  Figure  Cl. 

k]  forward  =  2.82e-09  exp  (-75.4  T  ^  )  (C5) 

k i reverse  =  2.82e-09  exp  (-75.4  T'1^  -  3083.7T'1)  (C6) 

Using  the  harmonic  ocsillator  approximation  for  vibrational  energy  level  spacing  in 
the  Schwartz,  Slawsky,  and  Herzfeld  (SSH)  (Ref.  C5)  treatment  of  vibrational  relaxation 
rates,  the  CO(2-l)  relaxation  rate  constant  is  estimated  to  be  twice  that  for  CO(  1-0).  No 
experimental  information  is  available  for  the  CO(2-0)  relaxation  rate,  but  it  can  be  assumed 
that  the  rate  constant  will  be  much  smaller  than  the  (2-1)  process  since  it  involves 
transferring  two  quanta  of  energy  per  collision.  The  coefficients  suggested  for  these  upper 
level  vibrational  transitions  are 

k 2 forward  =  5.64e-09  exp  (-75.4  T^ )  (C7) 

^reverse  =  5.64e-09  exp  (-75.4  T'1^  -3045.6T1)  (C8) 

k 3 forward  =  5.64e-10  exp  (-150.8  T>/3)  (C9) 
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^reverse  =  5.64e- 10  exp  (-150.8T'1/3  -  6129.3T’1) 


(CIO) 


The  appropriate  energy  level  spacing  has  been  used  for  the  temperature  defect  coefficients 
for  the  reverse  reactions. 
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APPENDIX  D:  V-T  ENERGY  TRANSFER  FROM  CO  TO  02  OR  N2 


The  following  vibrational-to-translational  energy  transfer  processes  are  possible  for  CO 
excited  in  the  v=l  or  v=2  levels,  when  CO  collides  with  an  atmospheric  modecule 


CO(l)  +  M  =  M  +  CO(O) 

(Dl) 

CO(2)  +  M  =  M  +  CO(l) 

(D2) 

CO(2)  +  M  =  M  +  CO(O) 

(D3) 

At  the  altitudes  discussed  in  this  review,  the  dominant  molecular  species,  M,  are  N2  or  02. 

Taylor  (Ref.  Dl)  reported  the  following  temperature  dependence  for  the  rate 
constant  of  the  reaction  shown  in  Eq.  (Dl): 

^forward  =  6.67e-8  exp(-208.3T‘1/3)  (D4) 

It  is  difficult  to  distinguish  the  V-T  contribution  to  the  effective  rate  of  relaxation  in 
a  mixture  such  as  CO/N2  or  C0/02,  where  competitive  V-V  processes  can  occur  to 
resonant  or  nearly  resonant  vibrational  energy  levels.  Experimental  data  are  often  only 
available  for  the  intermolecular  V-V  transfer  rates  since  these  are  generally  much  faster  than 
the  V-T  processes.  However,  the  latter  can  be  inferred  from  CO-atom  energy  transfer  rates 
and  from  CO/CO  V-T  experimental  data.  Taylor  used  three  primary  references  (Refs.  D2- 
D4),  which  provided  data  for  CO  relaxation  by  Ar,  and  CO  self-relaxation.  The  size  of  the 
error  bars  on  these  measurements  is  large  (e.g.,  43%  in  Ref.  Dl),  so  rates  inferred  by 
Taylor  (and  also  in  the  present  analysis)  are  roughly  order  of  magnitude  estimates.  In 
Figure  Dl,  the  reported  data  points  are  plotted  without  error  bars  for  the  purpose  of  clarity. 
More  recent  but  consistent  experimental  data  on  energy  transfer  from  CO  to  He  (Refs.  D4, 
D5)  are  also  included  in  Figure  Dl. 

The  temperature  dependence  function  of  the  CO  to  M  (M=  02  or  N2)  V-T  rate  will 
probably  lie  in  between  those  for  CO  relaxation  by  He  and  by  Ar,  based  on  the  relative 
masses  of  the  collision  partners.  From  Figure  Dl,  it  is  evident  that  the  CO  self-relaxation 
rate  also  falls  in  between  the  He  and  Ar  cases.  Therefore,  a  linear  Landau-Teller  function 
was  fitted  to  the  Kovacs  and  Mack  (Ref.  D3)  CO-CO  data  points  and  an  average  In  k 
intercept  from  the  other  data.  Die  final  temperature  dependence  is 

ki  forward  =  9.90e-09  exp(-168.1T-l/3)  (D5) 
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In  k  (k  in  cm3 moleculc^s*1) 


Figure  Dl.  V-T  energy  transfer  processes  from  CO(l).  CO/Ardata 
from  Millikan  and  White  (Ref.  D2),A.  CO/He  data  from  Millikan  and 
White,  (Ref.  D2)  A;  Maricq  et  al.  (Ref.  D7),0;  Schinke  and  Dierkson, 
(Ref.  D6),«;  Miller  and  Millikan  (Ref.  D5),E.  CO/CO  data  from 
Kovacs  and  Mack  (Ref.  D3),0. 
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Comparing  the  rate  constant  for  the  V-T  process  at  any  given  temperature  with  that  of  the 
V-V  CO  to  N2  or  O2  rate  constants,  it  is  evident  that  the  latter  are  considerably  faster  at  low 
temperatures  (e.g.,  200  K).  However,  as  the  temperature  increases,  the  difference  in  rate 
constants  for  the  two  types  of  processes  decreases.  For  example,  at  1500  K,  the  V-T 
relaxation  of  CO  (v=l)  by  O2  is  approximately  equal  to  the  V-V  relaxation  by  O2. 

Using  the  harmonic  oscillator  approximation  in  the  Schwartz,  Slawsky,  and 
Herzfeld  (SSH)  theorem  (Ref.  D8),  the  rate  constant  for  the  (2-1)  relaxation  of  CO  [Eq 
(D2)]  is  estimated  to  be 


k2  forward  =  1.98e-08  exp(-168.1TU3)  (D6) 

The  appropriate  CO  vibrational  energy  level  spacing  can  be  included  in  the  (-  AE/T)  term 
for  the  reverse  process.  The  (2-0)  relaxation  level  should  be  a  much  less  probable 
transition  since  it  requires  two  quanta  of  energy  to  be  transferred.  Therefore,  a  scaled 
estimate  for  the  (2-0)  rate  is 


^3  forward  =  1  -98e-09  exp(-336.2  T-^)  (D7) 

Table  1  of  the  main  text  contains  the  full  listing  of  rate  coefficients  for  Eqs  (D1  )-(D3), 
including  values  for  the  less  probable  reverse  processes. 
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APPENDIX  E:  CO  SPONTANEOUS  EMISSION  PROCESSES 


The  Einstein  A  coefficient  for  a  transition  from  an  upper  level  u  to  a  lower  level  1  is 
related  to  the  sum-of-lines  band  strength  Su|  by  the  following  expression  (Ref.  El): 

a  _  8tcv2Sui  exp[Ei/205.727]  ,Ei\ 

gu  l-exp[-v/205.727]  1  ' 

where  the  energy  of  the  standard  temperature  (296  K)  is  205.727  cm"1.  The  appropriate 
units  for  the  energy  of  the  lower  level,  Ei,  and  the  energy  of  the  transition,  v,  are  cm'1. 
The  units  for  the  speed  of  light,  c,  are  cm/s.  The  degeneracy  factor  of  the  upper  level,  gu, 
is  dimensionless. 

The  1987  Atmospheric  Trace  Molecule  Spectroscopy  (ATMOS)  sum-of-lines  band 
strengths  (Ref.  E2)  for  the  CO  vibrational  level  transitions  are  9.81  x  10*18  (1-0),  5.70  x 
10' 22  (2-1),  and  6.74  x  10*20  (2-0)  in  units  of  cm*1  /  molecule*1  cm*2  or  cm/molecule  (i.e., 
one  Benedict)  defined  at  296  K.  The  calculated  spontaneous  emission  rate  constants  are 
33.98  s*1,  64.43  s*1,  and  0.92  s*1,  respectively.  The  spacing  of  the  vibrational  levels  is 
given  in  Table  2  of  the  main  text. 

The  ATMOS  band  strengths  are  comparable  to  those  in  the  HITRAN  (Ref.  E3) 
listing  of  the  same  year  (1986).  The  HITRAN  band  strengths  are  9.81  x  10*18  (1-0),  5.69 
x  10*22  (2-1),  and  7.52  x  10*20  (2-0)  cm/molecule.  Several  papers  that  contain 
spectroscopic  measurements  and  calculations  pertaining  to  the  (1-0)  and  (2-0)  levels  of  CO 
are  listed  in  Refs.  E4-E7. 
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APPENDIX  F  (CONT.):  MODEL  ATMOSPHERE  USED  FOR  SIIARC  CALCULATIONS 
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